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ABSTRACT 

We compute the three-point temperature correlation function of the COBE 
Differential Microwave Radiometer (DMR) first-year sky maps to search for 
non-Gaussian temperature fiuctuations. The level of fiuctuations seen in 
the computed correlation function are too large to be attributable solely to 
instrument noise. However the fiuctuations are consistent with the level expected 
to result from a superposition of instrument noise and sky signal arising from a 
Gaussian power law model of initial fiuctuations, with a quadrupole normalized 
amplitude of f7 //K and a power law spectral index n = 1. We place limits on 
the amplitude of intrinsic three-point correlations with a variety of predicted 
functional forms. 

Subject headings: cosmic microwave background — cosmology: observations 
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1. Introduction 

NASA's COBE collaboration has reported the detection of anisotropics in the cosmic 
microwave background (CMB) radiation (Smoot et al. 1992; Bennett et al. 1992a; Wright 
et al. 1992 ). The amplitude of these anisotropics is substantially greater than the 
limits on possible systematic effects remaining in the data (Kogut et al. 1992), while 
the pattern of anisotropy is not correlated with any known Galactic or extragalactic 
foreground (Bennett et al. 1992a; Bennett et al. 1993). A useful characterization of 
CMB anisotropics is provided by the n-point correlation functions, the average product 
of temperatures evaluated at n points in the sky with fixed relative orientation. The 
COBE collaboration reported a positive detection of the two-point temperature correlation 
function (Smoot et al. 1992; Wright et al. 1992), and found it to be consistent with 
infiationary and other cosmologies. One of the basic predictions of infiationary models is 
that they produce a nearly scale-invariant power spectrum of primordial density fiuctuations 
with a Gaussian distribution and random phases. It follows that large angular scale 
CMB temperature anisotropics should have vanishing three-point correlations on average 
(though any particular cosmic observer will, in general, measure small non-zero three-point 
correlations due to cosmic variance). In contrast, there are competing theories of structure 
formation (eg., theories based on topological defects or late-time phase transitions) that 
are compatible with present CMB observations for which the distribution of density and 
temperature fiuctuations are non-Gaussian. Such models can give rise to non- vanishing 
three-point temperature correlations, thus a crucial test of the Gaussian nature of the 
primordial density fiuctuations is provided by this function. Smoot et al. (1994) discuss a 
number of additional tests that can be applied to the DMR sky maps to search for evidence 
of non-Gaussian fiuctuations. 

In this paper we compute the three-point correlation function of the first-year DMR 
sky maps and compare the results to detailed Monte Carlo simulations of Gaussian, 
scale-invariant power law models. We conclude that there is evidence for non- vanishing 
three-point correlations in the data, but at a level consistent with the uncertainty associated 
with cosmic variance. In the absence of a positive detection of intrinsic three-point 
correlations we can place upper limits on correlations with selected, theoretically motivated 
functional forms. Falk, Rangarajan, & Srednicki (1993) have considered the effect of a 
cubic term in the effective potential of an infiationary field and have concluded that 
the leading contribution to the three-point function can be expressed as a quadratic 
function of the corresponding two-point function. Luo & Schramm (1993) have argued that 
nonlinear evolution of the gravitational potential cf) will give rise to a similar form for the 
three-point function with an amplitude that likely dominates the infiationary contribution. 
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In addition, Luo & Schramm (1993) suggest that a late-time phase transition could produce 
a contribution that may be expressed as a cubic function of the corresponding two-point 
function. We place constraints on the amplitudes of such contributions. 

The DMR experiment (Smoot et al. 1990; Bennett et al. 1992b) has produced two 
independent microwave maps (A and B) at each of 3 frequencies (31.5, 53 and 90 GHz). The 
results presented here are based on the first year of data, and we consider only the relatively 
sensitive 53 and 90 GHz maps. In all cases we restrict our analysis to Galactic latitudes 
|b| > 20°, and we subtract a best-fit monopole and dipole from the data. The general 
three-point correlation function is the average product of three temperatures with a fixed 
relative separation on the sky: 6*3(^1,^2,^3) = {T {fii)T {fi2)T [fi^)) where hi ■ fi2 = cos^i, 
^2 ■ h^ = cos 62 and fi^ ■ hi = cos ^3. Due to computational constraints we evaluate only 
two special cases: 1) the direction vectors n^, n^, and hk form an equilateral triangle, 
01 = 02 = 03 (the equilateral case), and 2) hi is nearly parallel to hj (the pseudo-coUapsed 
case). We motivate and define the use of "nearly" in §3 below. Our notation is as follows: 
latin indices z, j, and k refer to sky pixels, while greek indices a and /3 refer to bins of 
angular separation, which we take to be of width 2.6 degrees, the width of our sky pixels. 
Additionally, the use of as an argument to a correlation function denotes a value of 
angular separation in degrees. 

2. The Equilateral Case 

The DMR sky maps are binned into 6144 nearly equal-area pixels, 4016 of which pass 
the 20 degree Galaxy cut. We define the equilateral three-point correlation function to be 

i,j,k 

where the sum is restricted to pixel triples (z, j, k) for which all three pixel separations reside 
in a single angular separation bin, Ti is the observed temperature in pixel i after monopole 
and dipole subtraction, and Wi is the statistical weight of pixel i: Wi = l/crf = Nij (Tq, where 
Ni is the number of observations of pixel z, and ctq is the rms temperature fiuctuation per 
observation due to instrument noise. 

We have evaluated this statistic for the 53 and 90 GHz DMR maps and present the 
results in Figure 1. The left-hand panels of the figure show results for both the (A-|-B)/2 
("sum") maps, which include both sky signal and instrument noise, and the (A— B)/2 
("difference") maps which contain only instrument noise. The grey band in these panels 
represents the rms scatter computed from 2000 Monte Carlo realizations of the instrument 
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noise. The right-hand panels repeat the results for the sum maps, but in these panels 
the grey bands represent the rms scatter computed from 2000 Monte Carlo realizations of 
simulated CMB sky signal plus instrument noise. In the present analysis we have simulated 
the CMB signal using a Gaussian, scale-invariant power law model with random phases and 
a mean quadrupole normalization of 17 fiK. Power from / = 2 to / = 40 filtered through the 
DMR window function (Wright et al. 1994) is included. Instrument noise was simulated for 
each channel A and B (with the appropriate DMR sky coverage pattern) using the following 
values for the instrument noise per 0.5 second observation: ctq = 24.2, 28.2, 47.8, and 35.6 
mK (thermodynamic temperature) for channel 53A, 53B, 90A, and 90B respectively. For 
each realization of sky signal and noise, we form sum and difference maps and subtract 
best-fit (for |b| > 20°) monopoles and dipoles from each. From the ensemble of realizations 
we compute the rms scatter in each of the 48 angular separations bins, as well as the full 
covariance matrix to account for bin-to-bin correlations induced by the sky signal. 

We have run a set of noiseless Monte Carlo simulations to compare the results of 
our routine with existing analytic predictions for the level of cosmic variance expected 
in the three-point function (Srednicki 1993). First we have generated a set of 400 
Harrison-Zel'dovich (n = 1) skies with unit quadrupole normalization including power from 
/ = 3 to / = 40, filtered through the DMR window function. For each realization in this 
ensemble we compute the three-point function using uniform weights over the entire sky. 
The rms scatter of these functions may then be directly compared to Figure 2 of Srednicki 
(1993); our results agree to within a few percent. The Galactic plane cut is an additional 
source of theoretical uncertainty (a form of "sample variance"). We find that restricting 
the sky coverage to |b| > 20° produces an rms that is larger than the full-sky case by a 
factor of ~ 1.5, somewhat larger than the factor of 1.23 predicted by Srednicki, based on 
the work of Scott, Srednicki, & White (1993). An additional source of sample variance 
in our data arises from the pattern of weights used to estimate the three-point function. 
Including the quadrupole in our noise-free simulations, and weighting the computation of 
the three-point function by the inverse of the DMR noise variance, we find the rms to be 
a factor of ~ 5 — 8 times larger (depending on the angular separation) than the uniformly 
weighted, full-sky, no-quadrupole case. It is advantageous to weight the computation of the 
three-point function to minimize the noise variance since the signal-to-noise ratio in our 
first-year maps is relatively poor. However, with additional years of data, cosmic variance 
will dominate noise uncertainties and it will be preferable to weight the pixels uniformly to 
minimize this source of uncertainty. 

The hypothesis that the computed three-point function is consistent with zero is tested 
by computing = C„ {M~^)^^j Cj3 where C„ is the observed three-point function in 
angular separation bin a, and Map is the covariance matrix computed from the simulated 
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three-point functions. The observed values of S'l'e given in Table 1, the numbers given 
in parentheses are the percentage of simulations for which exceeded the observed value. 
Consider first the hypothesis that the observed three-point functions are consistent with no 
sky signal, ie., that the observed fiuctuations are due only to instrument noise. In this case 
we evaluate x^ using the three-point functions observed in the sum maps and the covariance 
matrix derived from the ensemble of difference maps which have the same noise properties 
as the sum maps. We find x^ = 69 and 57 for the 53 and 90 GHz data respectively; only 
4% of our 53 GHz difference simulations and 18% of our 90 GHz difference simulations had 
higher values of x^ ■ We conclude that there exist non- vanishing equilateral three-point 
correlations in the observed CMB sky at the ~95% confidence level. We next test the 
hypothesis that the fiuctuations seen in the three-point function are consistent with the 
level expected to arise from a superposition of instrument noise and sky signal arising from 
Gaussian, scale-invariant power law fiuctuations. In this case x^ evaluated using the 
covariance matrix computed from the ensemble of sum maps; we find x^ = 49 and 50 for 
the 53 and 90 GHz data respectively. As indicated in Table I, these values are well within 
the range seen in the simulations, which indicates that the fiuctuations observed in the 
equilateral three-point function are consistent with a superposition of instrument noise and 
Gaussian CMB fiuctuations. 

As a check on the noise levels in our simulations we have evaluated x^ using the 
three-point functions observed in the difference maps. We find the 53 GHz value to be well 
within the expected range while the 90 GHz value is not: only I of our 2000 simulations 
had as large a x^ foi" 48 degrees of freedom. We attribute this to a weak anti-correlation 
between the noise in the 90 GHz A and B channels, which we have recently identified by 
other analyses. As a consequence the 90 GHz difference map is not a reliable estimator of 
the instrument noise level in the maps. We expect to correct this anti-correlation in future 
releases of the data and to discuss it more fully in a future publication. We emphasize that 
no present or previously reported results are compromised by this effect. 

In the absence of a positive detection of intrinsic three-point correlations we can place 
limits on the amplitudes of contributions with certain predicted functional forms. Falk et al. 
(1993) and Luo & Schramm (1993) argue that the equilateral configuration of the three-point 
function will have leading contributions of the form C^\9) ~ A;2C|(^) -|- A;3C|(^), where 
6*2 (^) is the corresponding two-point function, and k2 and are constants. To place limits 
on k2 and k^ one should ideally invoke a full Monte Carlo simulation of CMB skies with the 
appropriate non-Gaussian effects included. However, since the above predictions are only 
approximate, we feel it is sufficient at this time to adopt the theoretical form for the mean 
of 6*2 (^) as a basis for our model functions, and to perform the fits using the errors obtained 
from our Gaussian Monte Carlo simulations. In particular we define C2{0) as (Abbott & 
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Wise 1984) 

40 

C2(^) = (6/5)B2/ + 1)/W + 1)) G^W) 

1=2 

where the G/'s, given by Wright et aL (1994), define the DMR window function, and 
the P;'s are the Legendre polynomials. For reference, 6*2(0) = 4.24. We then perform 
least-squares fits to the data by minimizing a defined to be 

a,f3 

for m = 2, 3. We evaluate the model function C2{0) at the center of each angular separation 
bin. The covariance matrix is derived from the ensemble of sum maps; thus the quoted 
confidence intervals include the effects of cosmic variance. Since the two functional forms 
are nearly degenerate, we cannot place meaningful simultaneous constraints on the two 
parameters k2 and so we perform each fit separately. We define approximate 68% 
confidence intervals for k2 and ks by the condition x^(^m) < Xmin + 1 m = 2,3. The 
results of these fits are given in Table 2, which also includes weighted averages. Both k2 
and ks are consistent with zero. 

3. The Pseudo-Collapsed Case 

The simplest configuration of the three-point correlation function is that in which two 
legs of the correlation function are evaluated at a common point in the sky. SubbaRao et 
al. (1993) have proposed the use of the "collapsed" three-point function as an effective 
probe of non-Gaussian structure in the CMB. Using the same notation as above, we define 
the collapsed three-point function to be 

i,3 i,3 

where the sum is over all pixel-pairs within a given angular separation bin. Unfortunately, 
because of the quadratic term, Tj, in this statistic, the collapsed three-point function 

(c) 

suffers from a severe noise bias. In other words C3 is not a central estimator of the true 
collapsed three-point function in the presence of instrument noise. To see this consider 
an ensemble average of C^"^ over many realizations of the instrument noise, with a fixed 
true sky temperature. We can decompose the observed sky temperature Ti into a true 
temperature ti and a noise contribution n^: Ti = ti -\- rii. Then, by using the facts that 
{ti) = ti^ (rii) = (nf) = 0, and (nj) = af we obtain 

i,3 i,3 i,3 i,3 
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The first term is tlie "true" tliree-point function of tlie sky temperature, wliile tlie second 
is a bias term tliat is a cross-correlation of tlie true sky temperature with the average noise 
pattern. Given the current DMR noise levels and the apparent absence of three-point 
correlations in the CMB anisotropics, the bias term dominates the behavior of this statistic. 
We have investigated a number of ways around this problem. For example, we could select 
the weights Wi in such a way as to make the bias term cancel based on the average noise 
properties of the map. However, detailed Monte Carlo simulations show that no choice 
of weights simultaneously removes the noise bias while recovering the "true" three-point 
function. Alternatively we can evaluate the collapsed function with a subset of the maps, 
wherein we replace the quadratic term with the expression Ta,jTb,j where Ta is the 
observed temperature in the channel A map, and likewise for Tg. Since the noise in the 
two channels is uncorrelated (the 90 GHz channels notwithstanding, see above), there is no 
surviving bias. However one does sacrifice significant sensitivity with this scheme. Smoot et 
al. (1994) give a more detailed discussion of noise bias in higher-order statistical quantities. 

The technique adopted in this paper for overcoming the noise bias is to evaluate a 
"pseudo-coUapsed" three-point correlation function. We define this statistic as 

i,j,k 

where the sum on j is over all pixels that are nearest neighbors to z, and the sum on k is 
over all pixels (except j) within a given angular separation bin of i. Since the sky map 
pixels are 2.6° across and our beam has a 7° FWHM, the nearest neighbor pixels will have 
correlated sky signal but uncorrelated instrument noise. Additionally, since most pixels 
have eight nearest neighbors, we are summing over roughly four times as many independent 
pixel triples as in the true collapsed case, which significantly enhances the sensitivity. 

We have evaluated this statistic for the 53 and 90 GHz DMR maps and present the 
results in Figure 2, which has the same format as Figure 1. As with the equilateral case, 
we test the hypothesis that the data are consistent with vanishing three-point correlations 
by computing = lZa,fiCa {^~^)afj^P where C„ is now the observed pseudo-coUapsed 
three-point function in angular separation bin a, and Map is the covariance matrix 
computed from the ensemble of simulated, pseudo-coUapsed three-point functions. The 
values of are given in Table 3, which has the same format as Table 1. The values 
for the hypothesis that the three-point correlations in the sum maps are consistent with 
instrument noise are 138 and 127 for the 53 and 90 GHz data respectively. Only 0.1% of 
our 53 and 90 GHz difference simulations had a larger x^. We conclude that there exist 
non- vanishing pseudo-coUapsed three-point correlations in our sky at the ~99.9% confidence 
level. However, when we include the effects of cosmic variance, the values diminish 
considerably. Both the 53 and 90 GHz values fall well within the range of our simulations. 



- 9 - 



as indicated in Table 3. Thus cosmic variance can comfortably explain the observed level of 
pseudo-coUapsed three-point correlations. 

The functional form predicted by Luo & Schramm (1993) for the pseudo-coUapsed 
configuration of the three-point function is 

where we have used the fact that the mean separation between nearest neighbor pixels is 
~ 2.6°. For the form of C2 given in §2 we find 6*2(2. 6°) = 4.09. Using the above model 
functions and the covariance matrices computed from the ensemble of pseudo-coUapsed 
functions, we perform least squares fits for k2 and as in §2. The results are given in Table 
4. As before, the results are consistent with zero; additionally, since the pseudo-coUapsed 
configuration has greater sensitivity, the confidence intervals are somewhat smaller than 
those derived from the equilateral configuration. 

4. Conclusions 

We have evaluated the three-point temperature correlation function for the first year 
COBE DMR sky maps and find evidence for non-zero three-point correlations in the data. 
However, we demonstrate that the observed level of fiuctuations are consistent with the 
level expected to result from the superposition of instrument noise and CMB sky signal 
arising from a Gaussian, scale-invariant power law model of initial fiuctuations, with a 
quadrupole normalized amplitude of 17 //K. We have placed limits on the amplitudes of 
intrinsic three-point correlations with specific, theoretically-motivated functional forms. 

Given that the three-point function is a cubic statistic, the noise levels will diminish 
relatively rapidly (oc time~^^'^) with additional data. In the four year 53 GHz maps, the 
rms noise per 2.6° angular separation bin will be roughly eight times smaller than the 
levels depicted in Figures la and 2a. This corresponds to ~(11-12 /uK)"^ for the equilateral 
configuration and ~(6-9 /uK)"^ for the pseudo-coUapsed configuration, depending on angular 
separation. Thus, our sensitivity to three-point correlations will ultimately be limited by 
cosmic variance. 



Table 1. Values for the Equilateral Configuration' 



Data 



(-)^ 

af) I 



53 GHz C(+) 
53 GHz C(-) 
90 GHz C(+) 
90 GHz C(-) 



49 (42%) 

50 (39%) 



69 (4%) 
5f (36%) 
57 (f8%) 
87 (0.05%) 



^There are 48 angular separation bins in the equilateral three-point function. 
^C^a'^ denotes a three-point function computed from an (AibB)/2 map. 
■^Af^^^ denotes a covariance matrix computed from an ensemble of simulated (AibB)/2 
maps; separate covariance matrices were computed for the 53 and 90 GHz noise levels. 



Table 2. 68% Confidence Intervals, in //K"^, for and k-^ 



Frequency^ k2^ k-^^ 



53 GHz -324 <k2< 2841 -167 < < 750 

90 GHz -7325 <k2< -431 -2020 < < 118 

Average -1074 <k2< 1803 -324 <k^< 519 



^Approximate 68% confidence intervals, derived from the equilateral configuration, include 

uncertainties due to cosmic variance. 

^The line denoted "Average" gives weighted average confidence intervals under the 

assumption that the 53 and 90 GHz uncertainties are uncorrelated. 
^k^ is the coefficient of the model function C\{Q) 

^k-i is the coefficient of the model function C\{Q) 
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Table 3. Values lor the Pseudo-coUapsed Configuration^ 



Data'- X^(Mff) X^l^^)- 

53 GHz C(+) 50 (83%) 138 (0.01%) 
53 GHz C(-) • • • 71 (42%) 

90 GHz C(+) 89 (18%) 127 (0.01%) 

90 GHz C(-) ••• 114(0.9%) 



^There are 71 angular separation bins in the pseudo-coUapsed three-point lunction. 
^C^^^ denotes a three-point lunction computed Irom an (AibB)/2 map. 
'''M^j denotes a covariance matrix computed Irom an ensemble ol simulated (AibB)/2 
maps; separate covariance matrices were computed lor the 53 and 90 GHz noise levels. 



Table 4. 68% Confidence Intervals, in //K"^, lor k2 and 



Frequency^ k^'^ 

53 GHz -169 <k2 < QU -9 < A^s < 534 

90 GHz -1223 <k2 < 113 -1448 < < -353 

Average -323 < k2 < 371 -212 < ks < 276 



^Approximate 68% confidence intervals, derived Irom the pseudo-coUapsed configuration, 

include uncertainties due to cosmic variance. 

^The line denoted "Average" gives weighted average confidence intervals under the 

assumption that the 53 and 90 GHz uncertainties are uncorrelated. 
^k2 is the coeflicient ol the model lunction [26*2(2. 6°)C2(^) + C|(^)] 

'^ks is the coefficient of the model function 6*2(2. 6°)C|(^) 
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Fig. 1. — a) 53 GHz equilateral three-point correlation functions. The thick line is the result 
for the (A-|-B)/2 map, the thin line is for the (A— B)/2 map. The grey band represents the 
rms scatter due to instrument noise, see the text for details, b) 53 GHz equilateral three-point 
function for the (A-|-B)/2 map. The grey band represents the rms scatter due to simulated 
CMB sky signal plus instrument noise, see the text for details, c) 90 GHz equilateral three- 
point correlation functions, same format as a), d) 90 GHz equilateral three-point correlation 
function, same format as b). 
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Fig. 2. — a) 53 GHz pseudo-coUapsed three-point correlation lunctions, same format as 
Figure la. b) 53 GHz pseudo-coUapsed three-point lunction, same format as Figure lb. c) 
90 GHz pseudo-coUapsed three-point correlation functions, same format as Figure Ic. d) 90 
GHz pseudo-coUapsed three-point correlation function, same format as Figure Id. 



